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Abstract 

A phenomenological Dyson-Schwinger equation approach to QCD, formalised 
in terms of a QCD based model field theory, is used to calculate the electro- 
magnetic charge radius of the pion. The contributions from the quark core 
and pion loop, as defined in this approach, are identified and compared. It is 
shown explicitly that the divergence of the charge radius in the chiral limit 
is due to the pion loop and that, at the physical value of the pion mass, 
this loop contributes less than 15% to (r^); i.e. the quark core is the domi- 
nant determining characteristic for the pion. This suggests that quark based 
models which fail to reproduce the In 771,^ divergence of (r^) may nevertheless 
incorporate the dominant characteristic of the pion: its quark core. 
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Typeset using REVT^ 



I. INTRODUCTION 



Chiral symmetry and its dynamical breakdown is crucial to the nature of the strong 
interaction spectrum and the study of this in Quantum Chromodynamics (QCD), and models 
thereof, is of great current interest. The most striking signals of Dynamical Chiral SjTiimetry 
Breaking (DCSB) are found in the kinematic and dynamical properties of the pion and one 
tool that has been used extensively to study these properties is Chiral Perturbation Theory 
(ChPT) Chiral Perturbation Theory can be understood as the study of the necessary 

consequences of the chiral Ward Identities via the construction of an effective action, using 
field variables with pionic quantum numbers, in such a way as to ensure that these identities 
are realised. It should be noted that in this approach the pionic field has no physical 
significance - it is merely an auxiliary field - and should not be identified with the physical 
pion 0. 

At first nonleading order in ChPT, O(E^), the effective action is only completely de- 
termined once the effect of one-pion loops, generated by the O(E^) part of the action, is 
included. The regularisation of the divergence of each of these loops introduces ten arbitrary 
parameters at this level. The theory is not renormalisable but the approach is nevertheless 
useful because higher order loops do not contribute to the O(E^) action 0; i.e., the action 
at O(E^) is completely specified once the effect of the O(E^) loops is taken into account and 
the 10 parameters fixed by comparison with experimental data. 

The pion (or, more generally, pseudoscalar) loops in ChPT are characteristic of the 
approach and, indeed, are sometimes regarded as being the dominant feature. The expression 
for every physical observable receives a contribution from such loops which depends on the 
mass of the particle in the loop and which diverges in the chiral limit. For example, in the 
case of the electromagnetic pion charge radius the loop introduces a logarithmic dependence 
and hence a logarithmic divergence in the chiral limit. Herein we are interested in estimating 
the importance of these loop contributions, evaluated at real pseudoscalar masses, relative 
to what one might call the "quark core" contribution and our exemplary case will be the 
electromagnetic charge radius of the pion. 

The pseudoscalar mesons are the simplest bound states to study within the coupled 
Dyson-Schwinger-Bethe-Salpeter equation approach to QCD, recent reviews of which can 
be found in Refs. In vector exchange theories with DCSB, such as QCD, the Goldstone- 

boson character of the pseudoscalar excitations can be understood through the existence of 
an identity between the fermion Dyson- Schwinger equation (DSE) and the pseudoscalar 
Bethe-Salpeter equation (BSE). If, in the limit —>■ 0, this BSE reduces to the DSE then 
DCSB necessarily entails the existence of massless pseudoscalar excitations. This result was 
established for the rainbow-DSE-ladder-BSE pair in Ref. [^]. (The rainbow approximation 
in the DSE is equivalent to the ladder approximation in the BSE.) In this approach, the 
"quark core" of a meson is described by the Bethe-Salpeter amplitude obtained as a solution 
of the pesudoscalar BSE; the kernel of which, in practice, is often constructed as a generalised 
ladder kernel: generalised in the sense that the quark propagators are obtained as solutions 
of the rainbow DSE with a gluon propagator whose large spacelike-g^ behaviour is that given 
by the renormalisation group in QCD and whose small spacelike-g^ behaviour is modelled 
so as to ensure confinement [p!T|-p!^. 

The coupled rainbow-DSE-ladder-BSE approach to modelling QCD is equivalent to the 
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Global Colour-symmetry Model (GCM), first described in Ref. [|T5|] and extended to the 
pseudoscalar sector in Refs. |lT6ip!7| . The extension to other mesons and baryons is reviewed 
in Ref. |T^. As discussed in Ref. Jl^, this approach provides a very good description of 
kinematic and dynamical properties of pions. For example, using the quark core alone, it 
describes tt-tt partial wave amplitudes very well up to 600 MeV; i.e., from threshold up to the 
first resonance. As indicated above, the parameters in the GCM describe the quark-quark 
interaction at small spacelike-g^. 

The model represents QCD as an effective two-point theory, coloured quark currents 
interact via dressed gluon exchange, and has the "disadvantage" that there is some model 
dependence in the results: for example, it predicts values for the parameters that appear in 
the effective lagrangian of ChPT |ll7| , [T9| . However, since these parameters have been fitted 
to low energy phenomena in QCD, this provides a means of checking the model assumptions 
in the GCM. Thus the model dependence is also an advantage: the assumptions on which 
the model is founded can be tested and this enables the development of an understanding 
of hadronic phenomena in terms of quark and gluon interactions. 

Herein, in the context of the GCM, we calculate the electromagnetic charge radius of the 
pion taking into account both the "quark core" and the pion-loop rescattering correction. 



The "quark core" contribution, obtained in a generalised impulse approximation [0,ED|,BT 



is regular in the chiral limit whereas, as noted above, the pion-loop contribution is not. Our 
aim is to estimate which is the dominant piece at the physical pion mass and hence to infer 
whether the electromagnetic structure of the pion is dominated by the "quark core" or a 
"pseudoscalar meson-cloud" . 

In Sec. II we describe the calculation of {r"^)^^^, the generalised impulse approximation 
"quark core" contribution to the charge radius. In Sec. HI we discuss the calculation of the 
pion loop contribution to the charge radius. We put these pieces together to obtain the net 
7C7['~f vertex and obtain an expression for the net charge form factor in Sec. IV. Our results 
are presented and discussed in Sec. V. We summarise and conclude in Sec. VI. 



II. GENERALISED IMPULSE APPROXIMATION 



In Euclidean space, with metric 6^^, = diag(l, 1, 1, 1), 7^ hermitian and = m^, we 
define the generalised impulse approximation contribution to the connected inT'-f vertex to 



be: 17.20.2^ 



(1) 



J 77 



(2^)4 



T^{ko+]p)S{k++)iTf,{k++, A;__)5'(A;__)r^(A;_o; q)S{k^^ 



In this expression the trace over colour and flavour indices has been carried out leaving only 
the trace over Dirac indices and the definition 



k. 



af3 



k + ^p+^q 



(2) 



has been used. The dressed quark-photon vertex is denoted by T^{pi,p2), the pion Bethe- 
Salpeter amplitude by r^(p;P) and the dressed quark propagator by S{p). This vertex is 
illustrated in Fig. ^ and is what we mean by the "quark core" contribution to the electro- 
magnetic properties of the pion. 
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A. Quark propagator 



We describe equation (|ip as "generalised impulse approximation" for a number of reasons. 
First, the internal quark lines are described by dressed quark propagators that are obtained 
by solving the quark Dyson- Schwinger equation (DSE): 



S-\p) =r^.p + m+y^ j ^^^S{k)n{k,p)D,,{{p - kf). (3) 

In this equation D^u{k) is the dressed gluon propagator and r^(pi,p2) is the dressed quark- 
gluon vertex; each of these quantities satisfies its own DSE. The general solution of this 
equation has the form 

Sip) = -11 -pavip^) + cTsip^) . (4) 

A commonly used equivalent representation is to specify the inverse in terms of two scalar 
functions: 

S-\p) =i-f ■pA{p'^) +m + B{p'^) (5) 

where m is the current quark mass. Both forms will be used herein. The commonly used 
"rainbow approximation" DSE is obtained from Eq. by writing 

nik,p) = iu. (6) 

Equation (^) has been studied extensively using many model forms for the dressed gluon 
propagator and quark-gluon vertex and some of this work is summarised in Refs. A 
good deal has been learnt about the analytic properties of the quark propagator in the 



complex plane [^,0 and this is important in the calculation of many observables. In 
connection with the electromagnetic charge radius of the pion, however, it is sufficient to 
have a representation (or interpolation) of the quark propagator that is valid on the real 
spacelike axis. One such simple form, which is a modification of that used in Ref. and 
a simplification of the form used in Ref. |T^, is 



Mp') = 1 , (7) 
p' + Al 



B{P') = zri^ , (8) 



where reasonable values of the parameters are Ai ~ 0.5 GeV ~ A2. We remark that the fact 
that B 0, even when the quark bare mass is zero, is a manifestation of DCSB and, in 
general, for the small quark bare masses that are appropriate for the pion sector, ^(p^) and 
Bip^) are insensitive to m; a result that can be seen clearly in the model of Ref. [^. This 
means that in our parametrisation of the DSE results Ai and A2 are also insensitive to m. 
We emphasise that these forms are only suitable for interpolation on the real spacelike axis. 
We will discuss this Ansatz in more detail below [see Sec. |1I C|| . 
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B. Pion Bethe-Salpeter amplitude 



The generalised impulse approximation, in addition to using dressed quark propagators, 
uses a pion Bethe-Salpeter amplitude that is obtained as the solution of the homogeneous 
Bethe-Salpeter equation (BSE): 

K'iP;P) = J ^K^'-''''{k,p;P) (5(A;-iP)r.(A:;P)5(A: + iP))*" (9) 

where P is the centre-of-mass momentum of the bound state, p is the relative momentum 
between the quarks in the bound state, and the superscripts, which run from 1-4, are asso- 
ciated with the Dirac structure of the amplitude. The kernel, K^'^'^^{p, k; P), has no flavour 
structure in the isospin symmetric limit = and since T.„{p; P) oc Ic and S{p) oc Ic 
then K^^'^'^{p, k; P) oc Ic, where Ic is the identity in colour space. The ladder approximation 
BSE is obtained from Eq. (P) by writing 

P) = lg'D,,{p - k) {^,r ■ (10) 

Here D^i, is a model dressed gluon propagator which, in Landau gauge, has the form: 

^D,Ak)=(^6,.-'^^D{k^). (11) 

The behaviour of D{k^) at large spacelike- A;^ is determined from the QCD renormalisation 
group: 



Avr 

k^ In 



Dik') ^ -.^^ (12) 





^^QCD 

with A = 12/[33 — 2Nf] and where Aqcd ~ 0.2 GeV is the renormalisation point invariant 
QCD mass scale. The behaviour of D{k^) at small and intermediate spacelike-fc^ is unknown 
and in the GCM [or coupled DSE-BSE approach] a few parameters are used to specify a 
model form; i.e., the behaviour of the quark-quark interaction in the infrared is parametrised. 
The results one obtains for physical quantities in this approach are completely determined 
by the values of these parameters, as emphasised in Ref. Jl^. The constraints and commonly 
used parametrisations are discussed in detail in Refs. and one efficacious form has 

D{q') = CAlcB^'iQ) + , X (13) 

g2 In 

with C = (37r)3 and r = 3 fTT 



The most general form of the pseudoscalar amplitude, allowed by Lorentz covariance, 
which is odd under parity transformations is ||2^ 



T^ip; P) = 75 { E{p, P) - z(7 ■ p) (p ■ P) P(p, P) - z(7 ■ P) G(p, P) - [7 ■ p, 7 ■ P] H{p, P)} 

(14) 
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and since tt'' is even under charge-conjugation then E, F, G and H are even functions of 
{p ■ P) in this case. 

The coupled rainbow-DSE-ladder-BSE system, Eq. (|) with (|^) plus Eq. with (|10|), 
which is equivalent to the GCM effective action approach to QCD, has been studied ex- 
tensively and provides an excellent description of light-light, heavy-light and heavy- heavy 
meson systems [pJ] -[l4|. These studies show that it is a good approximation to write 



(15) 



In the chiral limit; i.e., when the current quark mass, m, is zero, the pseudoscalar 
ladder-BSE and quark rainbow-DSE are identical |10| and one has a massless excitation in 
the pseudoscalar channel with 



E(p,p2 = 0) oc 



(16) 



if the quark DSE allows a dynamically generated mass function: -B(p^) 7^ 0; i.e, if there is 
DCSB. This is the manner in which Goldstone's theorem is realised in the GCM and coupled 
DSE-BSE approach. In this case the parametrisation of Eq. (^ is also a good representation 
of the pion Bethe-Salpeter amplitude. 

The normalisation of the amplitude is fixed, as usual, by requiring that, for = — m^: 



(27r)4 {2ttY 



tr 



D 



+ 



(17) 



Nr. 



(2^)4 



tr 



D 



T^{k- P)S{h.)T^{k- -P)^^^ + r.(A;; P)^^^T^{k- -P)S{k,, 



dPu 



dPa 



If the Bethe-Salpeter interaction kernel, K{k,p; P), is independent of the centre-of-mass 
momentum, P, which is the case in ladder approximation, for example, then this reduces to 



Nr 



(18) 



(2^)4 



trz) 



K{k; P)Siko.)K{k; _p) ^^(^o+) + Y^{k- P) ^ V .(fc; -P)S{ko, 



dp, 



dp, 



C. Quark-photon vertex 

The generalised impulse approximation contribution also involves the dressed quark- 
photon vertex r^(pi,p2) which satisfies its own integral equation. Solving this equation is 



a difficult problem which has only recently begun to be addressed [^]. However, it is clear 
that if the quark propagator has momentum dependent dressing then r^(pi,p2) cannot be 
the bare vertex since the Ward Identity 

iPi'P2),tT,{puP2) = S-\pi) - S-\p2) (19) 

would not then be satisfied and hence the electromagnetic current of the pion would not be 
conserved. 
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It is possible to learn much about this vertex without actually solving its integral equa- 
tion. Indeed, it is argued in Refs. [^7|-30| that a "physically reasonable" Ansatz for this 
vertex must take the form: 



where |^| 



rf(p,A:) 



and 



+ 



r>,fc) = rf (p,fc) + r:(p,fc) 



(20) 



7m 



p2 _ 



[7 ■ p + 7 ■ fc] 



(21) 



{p-k)f,Tl{p, k) 



with rj(p,p) 



0. 



(22) 



In the bare quark propagator limit, A = 1 and B = constant, = 0. This Ansatz is 
"physically reasonable" because it satisfies the criteria a)-d) discussed in Ref. |Q: a) it 
satisfies the Ward-Takahashi identity; b) it is free of kinematic singularities; c) it reduces 
to the bare vertex in the free field limit; d) it has the same transformation properties as 
the bare vertex under charge conjugation and Lorentz transformations. The r^'^[p, k) piece 
is completely determined by the quark propagator but the remaining piece, Tj^{k,p), is 
undetermined. 

It is in this additional transverse piece that photon-vector-meson mixing contributions 
will appear; i.e., vector meson dominance contributions are confined to r^(p. A;). This is 

BC/ 



(23) 



especially clear in the Nambu-Jona-Lasinio (NJL) model where {p, k) = 7^ |^ 
We will define our generalised impulse approximation as Eq. (P with 

Tlip,k)^0. 



This allows for a very good description of the electromagnetic form factor in the spacelike 
region; i.e., away from resonance contributions |2l|. We note that criterion c) is important 
because, in association with a solution of the DSE for the quark propagator which exhibits 
the correct perturbative QCD leading- log behaviour of S(p), which is guaranteed when the 
gluon propagator in Eq. (^) has the ultraviolet form given in Eq. (p!^), it ensures that the 
vertex Ansatz will yield the same large-Q^ behaviour for the pion form factor as is predicted 
by QCD 0. 

With a quark-photon vertex that satisfies the above constraints and using the identities: 



— T 

r 



S{-kf = C^ Sik)C , 
-k;-p) = -C^T^{k;p)C 



T'^{-k,-p) = -C' T^{p,k)C, 
where C = 7274 is the charge conjugation matrix, it is easy to establish that 

{p + q),A^'\p,q)=0; 



(24) 
(25) 
(26) 



(27) 
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i.e., the 7r-current is conserved. Hence one can write 



GIA/ 



where 



Pp^q -(p + q)- QpiP -(p + q) 



(p + 

In the case of elastic scattering = and therefore Eq. (pSf ) reduces to 



Using Eqs. (p4| - p^ ) one can also show that A^^^(P, —P) is 



(28) 



(29) 



(30) 



2P,F 



GIA 



P 



(2^)4 



tr 



D 



d 



T^{k;P)S{k^o)T^{k;-P)—-S{k 



dp, 



d 



+T^{k; P)^S{k_o)T^{k; -P)S{k+o) 



(31) 



Comparing this with Eqs. (0) and (0) one observes that in generahsed impulse approxima- 
tion F'-'^^(0,P^) = 1 only if the Bethe-Salpeter kernel is independent of P; i.e., generalised 
impulse approximation combined with a P-independent Bethe-Salpeter kernel provides a 



consistent approximation scheme. In this case one has, in the chiral limit P 



-mz 



J TT 

Nr. 



ds s T^{sY (o-y - 2 [(Tsa's + sayCTy] 



-- 0: 
(32) 



crycry 



This illustrates that Eq. (|1]) is regular in the chiral limit and only weakly dependent on m^r. 

We are now in a position to consider the Ansatz for the quark propagator, Eqs. (0) and 
(§), in more detail. Following Ref. |]T^, this Ansatz for the fermion propagator can be used 
to fit the characteristic parameters of the vr-vr sector with the result that, with 



Ai = 0.527 GeV and A2 = 0.573 GeV, 

one obtains 

U = 0.093 GeV (0.093), = 0.53 fm (0.66) 
and the following values for the dimensionless scattering lengths: 



a° = 0.16 (0.20), 
a\ = 0.029 (0.038) 



-0.047 (-0.037), 
0.0014 (0.0017), 



-0.00037 



(33) 



(34) 



(35) 



where the experimental values, when known, are given in parentheses. These results were 
obtained using the approximation of Eq. (|T5[) and was calculated using Eq. (^2f). The 
formulae for the remaining quantites are all given in Ref. [17| and, apart from (r. 



2\GIA 



m 



Eq. ( |A1| ), will not be reproduced [^, however, all of the information needed to derive them 
is contained herein. The results are simply presented so as to demonstrate the efficacy and 
utility of this simple Ansatz and to fix reasonable values of the parameters. 
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D. Approximate form of F^^^{p'^,p- q,q'^ 



All of the elements of the calculation of F'^^^ijP'^p ■ q, q^) have now been described. The 
actual calculation involves a three-dimensional integral that can be evaluated numerically 
without difficulty for arbitrarily large spacelike momentum transfer provided a sensible DSE- 
BSE framework is used to describe the quark-gluon substructure of the pion. This is the 
subject of Ref. (The fact that Eqs. (^) and (H) should only be used for interpola- 

tion is important in this context.) Herein we are interested in the charge radius as it is 
affected by this parton substructure and by meson self-dressing. The detailed structure of 
F^^^ijP'^p ■ q, q^) is not necessary for this calculation. 

In general F^^^{p'^,p- q,q'^) depends on all three of its arguments and, as will become 
clear below, this dependence will be sampled by the pion loop contribution to the charge 
radius. However, our calculation is greatly simplified if we make the approximation that 



F^'\p\p.q,q'] 



(36) 



where = {p + q)"^, for then this piece appears only as a multiplicative factor in the 
pion loop contribution to the charge radius. This approximation is expected to lead to an 
overestimate of the importance of the pion loop to the charge radius since any off-shell effects 
in F^^^{p'^,p ■ q, q^) will lead to additional damping in the pion-loop integral, which is over 
spacelike momenta. 

Subject to this approximation then, for our purposes, all we need to know is 



d 



(37) 



Q2=o 



which has been presented in Ref. |]T^ and which we reproduce in Eq. ( [AID of the appendix. 
As we have seen, the simple interpolating form for the quark propagator that we are us- 
ing herein yields the charge radius in Eq. ( |5^ ) while the more sophisticated analysis of 
Ref. [0, using a somewhat different form of the fermion propagator, obtained the result 
= 0.59 fm. Both of these calculations are to be compared with the experimental value of 
r,r = 0.66 fm [^. This illustrates what we believe to be a general result: that the generalised 
impulse approximation consistently understimates the charge radius by < 15% leaving room 
for only small contributions from meson-dressing effects such as the pion loop and p'^-photon 
mixing. 



III. PION LOOP CONTRIBUTION 

We approximate the pion loop contribution to the connected 7r7r7 vertex as 

A;^L-P(p,g)= (38) 
-\ e'"' / WF ^''''(^ + q,p-k,k-q) DM DAp + q - k) A';'\k,p + q - k) 

where D^^i^k) is the pion propagator, A^^^{p,q) is given in Eq. (|ip and the vr-vr scattering 
kernel is T'^-'''''{pi + P2,Pi + P3,Pi + Pa), which can be written in terms of a single scalar 
function: 
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T^^'^Pl + P2, Pi + P3, Pi + Pi) 



(39) 



5'^^^ A{pi,p2; P3, Pa) + S'^'d'' A{pi,ps; pa, Pa) + 5'' 5'^ p^- ps, p^) 



with 



^(Pi,P2;P3,P4) = A{p2,Pi;p3,Pa) = A{pi,p2]Pi,pz) = A{p2,pi;p4,P3) (40) 
from which all of the relations of crossing symmetry follow. We also note that the result 

{p + q),A]f°np,q) = (41) 
follows from Eq. (|27D and as a consequence one may write 

A;^'^°°p(p, q) = T,{p, q) F^°°^{p\p- q, q') (42) 
which, in the case of elastic scattering [p^ = g^], reduces to 

A^^^^^lP, q) = e''' {p - q),F^'""'{iP + • (43) 

Equation (|38|) is illustrated in Fig. ||. This is the analogue in the GCM [and coupled DSE- 
BSE approach] of the one pion loop diagram calculated in ChPT. 



A. TT-vr scattering amplitude 



The TT-TT scattering amplitude is an important part of Eq. (E 



The form of 



A{pi,P2',P3,P4) near threshold in the GCM was studied in detail in Ref. fT^- Therein it 
was shown to reproduce the Weinberg term |^ at O(E^) and to possess additional struc- 
ture at O(E^) that provides for a better description of the scattering lengths and partial 
wave amplitudes in vr-vr scattering: 



A{s, t, u) 



ml + 2s — t — u 



(44) 



4iVc 



Ki (-I2ml + Qml{s + t + u) + 2 - - - 2{st + su+ tu 
+K2 [-2ml + s) {-2ml + t + u) + K3 (^-2ml + ml{s + t + u) - t 



where s = —{pi +^2)^, t = —{pi + ps)'^ and u = —{pi + Pa)^ are the usual Mandelstam 
variables and in the GCM 0: Ki = 0.000508, K2 = 0.00930 and = 0.00101. 

In the present study, however, an expansion in powers of momenta about threshold is 
inadequate since the integration in Eq. (^) samples momenta well away from threshold. 
A broad ranging discussion of the constraints that analyticity and unitarity place on the 



asymptotic behaviour of the scattering amplitude can be found in Refs. [37]; one result being 
that the forward scattering amplitude, A{s, 0,4m^ — s), cannot fall faster than 

The dominant contribution to A{pi,p2]P3,P4) in the GCM is obtained from Eqs. (A. 14) 



and (A. 16b) in Ref. ITi: 



d^^d'y ^tr [d.Uix) d.Uiy)^] fix - y) 



(45) 
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where, as usual, U{x) = exp(z7r ■ t/Ztt) and 

(2^) 



/>)= l&^^'''^fik') (46) 



with 

flp'f{p) = 2N,j ^,T^{efK{k+)K{k^) [k^Al + k.pA^A^ + \p'Al) , (47) 
[k± = k± [A± = A{k+) ± A{k-)] and where 

^(p) = p2A{p^y + 5(p2)2 • 

[This expression has been generahsed here to allow for momentum dependence of A[p^).] 
In Eq. (^71), the right-hand-side contributes the leading term in the expression for in the 
chiral limit: the other terms that contribute to Eq. ( ^2]) are obtained from Eqs. (A. 16a) 
and (A. 16c) in Ref. and in the point-meson [or Nambu-Jona-Lasinio model] limit of the 
GCM, defined such that Alp^) = 1 and B{p'^) = constant, these terms vanish. Taking this 
into account we renormalise Eq. (^) and define 

p' f{p) = }f,j^2 K{k^) ^(^-) (^'^- + ^ ■ + ¥ ^1) (49) 

with Af' chosen so that /(p^ = 0) = 1. 

The contribution of Eq. (H3) to v4(pi,p2;P3,P4) at tree level is easily found to be 



^'^°(Pl,P2;P3,P4) = 3^ (^3 (pi +p2)V(pi +P2) - EP?/(p')j • (50) 

We refer to this expression as the "leading-order" scattering amplitude because it provides 
the leading contribution at large s; i.e., it is the contribution to the scattering amplitude in 
the GCM which falls least rapidly with increasing s. We note that neglecting the momentum 
dependence of /(p); i.e., at lowest order in a momentum expansion, Eq. (|50D yields, upon 
continuation to Minkowski space, 

ALO(pi,P2;P3,P4) = ^^^1^ (51) 



which reproduces the result of Ref. |]35| in the chiral limit, — > 0. 

We would like an accurate estimate of this amplitude in the GCM. It is not difficult 
to evaluate the integral numerically using the results of DSE studies, as represented by 
Eqs. (1^, (^ and (|33|), but it is more instructive to distill these ingredients and obtain an 
analytic expression that manifests the important features. To do this we use Eqs. (^, (|^) and 
(p!5|) but, in addition, we make a further simplifying assumption, based on the "constituent 



mass" concept described in Ref. ||38|. The basic observation in this connection is that the 
integration in Eq. ( p7D is weighted, and cut off, by the factor k'^V.„{k'^Y that arises because 
of the quark-gluon substructure of the pion in the GCM. This manifestation of the quark 
core of the meson is a general feature of the coupled DSE-BSE approach and ensures that 
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even meson loop integrals are intrinsically finite. This factor concentrates the integration 
domain around s = A^, using Eq. ([T5|), and hence it is a reasonable approximation to replace 
B{p'^) in K{p^) by an "effective constituent quark mass": 



p2 + M2 

where a first estimate of M^. is obtained by evaluating B{p'^) 
mum of k"^ VT.ihP')'^: 



(52) 



at 



A2, which is the maxi- 



AL 
2A2 ■ 



(53) 



Using the values of Ai and A2 in Eq. (^ this yields ~ 0.22 GeV. 

A more sophisticated estimate of Mc may be obtained by requiring that using Eq. (|52| 
in Eq. (^) should yield the same result for /^r as the exact calculation. Using the approxi 
mations of Eqs.(|^, (H), ( [T5| ) one can rewrite Eq. ( |32D as 

N, 



fl = ^,l dss K{sf B{s) {b{s) - \B\s 



(54) 



which, not unexpectedly, is the result presented in Ref. Substituting Eq. (|52| ) for K in 
this expression one obtains 



J TT 



AT A6 (-7 A^ + 20 Ai M2 - 13 + In ^ [-4 A| + 4 A^ + 6 M- 

C 1 \ I\-2 



167r2 (M2 - Ai)4 

and requiring that this yield = 0.093 GeV gives 

M, = 0.31 GeV , 

which shows that Eq. (|53|) provides a reasonable estimate. 
Using Eqs. (0), (||), (|15|) and (|^) we obtain 



(55) 



(56) 



2 r°° 



;^ + Ai)2 



(57) 



with 



{M'c+^D ln^-2(M2-A2) 



(M2 - A2)3 

It is not difficult to establish that, at large-p2^ 



1 8 , p2 

— ;ii^T2 



(58) 



(59) 



Np^ A'i 

which gives A^'~'{s,t,u) consistent with the bounds of Refs. [0. We emphasise that A2 is 
not an arbitrary parameter introduced to regularise the integral but, rather, it is a direct 
measure of the quark-gluon substructure in the model arising as it does through solving the 
Bethe-Salpeter equation; i.e., it is a calculated quantity in the GCM. We note also that Ai 
has not disappeared but is contained implicitly in M^. 
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IV. NET 7r7r7 VERTEX 



The net niT'j vertex is obtained by adding Eq. directly to Eq. ([I|): 

A^ip, q) = e''^^Kf\p, q) + A;^-^°°p(p, q) (60) 
= e3*^T^(p,g)F^(p2,p-g,g2) (61) 



with 



h{p\p- q, q') = F^'\p',p- q, q') + F^^^^P^P' • (62) 



This does not lead to double-counting because the pion final-state-interactions cannot be re- 
cast into an additional dressing of the quark-photon vertex; unlike the p-meson contribution 
that is contained in T^{p, k) and which we neglected. Adding Eq. ( |5^ ) is thus a correction 
beyond generalised impulse approximation and as such leads to a renormalisation of the 
pion field in our approach, which is simply a reordering of the summation of diagrams in a 
well defined and systematic manner as prescribed by the bosonisation of the GCM. 

To clarify this statement we remark that going beyond generalised impulse approximation 
in the calculation of the electromagnetic pion form factor builds additional structure into 
the pion. In the framework of the GCM, the "bare pion"; i.e., the one whose Bethe-Salpeter 
amplitude appears in Eq. (|T|), is a ladder q-q bound state: it represents the "quark core" 



of the pion. As shown in Ref. [|17|, this core provides a very good description of kinematic 
and dynamical properties of the pion away from resonance contributions. Adding the pion 
loop contribution is the first step in allowing this ladder q-q bound state to dress itself 
with a cloud of ladder mesonic bound states. In the framework of the BSE, Eq. (^, this 
corresponds to going beyond ladder approximation. Indeed, this modification could be built 
into the pion BSE and would lead to a dependence on the centre-of-mass momentum, P, in 
the kernel of Eq. (^ thus changing the normalisation obtained from Eq. (|1^. This illustrates 
the nature of the GCM approach. The "tree-level" fields in the effective action represent 
ladder q-q bound states and the remaining interaction terms provide for meson self-dressing, 
in addition to meson-meson interactions, which extends the nature of the structure of the 
bound state beyond that of ladder approximation. 
The pion propagator in Eq. (BHI) has the form 



where E^(/c^) is the pion self energy which can be obtained as the solution of a DSE or using 
the eigen- value procedure of Ref. [|TS|. Herein we follow Ref. and use the approximation 



where (m^)^ — S^(— (m^)^) = 0. We have labelled the mass of the pion in the loop in 
order to artificially distinguish it from the mass of the external pions; a convention we adopt 
simply so as to make explicit the dependence of the charge radius on the mass of the loop 
particle. Equation (Bl) is expected to be a good approximation because the Bethe-Salpeter 
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amplitudes in Eqs. (P and (|38D ensure that the dominant integration domain is at small 
k'^ where the pion propagator should not be much modified from its on-shell form. Indeed, 
the dominant dressing contributions for the virtual pion propagator would come from a Svr 



intermediate state and this was shown to be negligible in Ref. 
Using Eq 



( |6^ in Eq. (^) we obtain the following general expression for the one-pion- 
loop corrected approximation to the electromagnetic pion form factor: 



1 f d'^k \ q ■ {p + q) k ■ {p + q) — q ■ k {p + 



(65) 



(2vr)^ 



p2g2 



{p ■ 



X 



A{p, -k;q,k-p~ q) - A{p, k-p-q;q, -k) pciAfp . 
(P + (mL)2)((-A;+p + g)2 + (r„L)2) - ^ ' 



-k +P + q), {-k + P + qf 



As we have mentioned above and make explicit below, the pion loop contribution is given 
by a finite integral because the "bare meson" in the loop already has an internal structure: 
its q-q core as given by the BSE, Sec. |11B. 



V. PION CHARGE RADIUS 



Equation ( p5[) is simplified by the approximations we have described above. Using 
Eq. (p6D, the "leading order" scattering amplitude, Eqs. (|50|) and (|57D, shifting the inte- 
gration variable: [/c — >• /c + g] , writing 



P = p — q and Q = p + q 



(66) 



introducing hyperspherical polar coordinates and placing the external pions on the chiral- 
limit mass-shell, [P^ + Q"^ = 0], we obtain the following expression ||0|: 



FAQ') = F^'\Q') { 1 + 



167r3/2 



3 



dk 



^dy (l + 2.,A 



(67) 



X 



(P - 2iQky - Q^)f{P - 2iQky - Q^) - ef{k^ 



[m. 



L\2] 



X In 



kQ{k-^-ikQy 
e - ikQy + (mL)2 + kQ^l^ 



fc2 — ikQy + (m^)2 — kQy/1 — y'^ 



Equation ( |67D can be used to calculate the one-pion-loop corrected value of the pion 
decay constant, which we will denote f^,, 



fl = flFM- 



A little algebra yields 

where, with f{x) given in Eq. ([ 



du 



[Au + {m^ 



L\2]2 



9{u) 



(68) 



(69) 
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g{u) = fiAu)+u-^f{Au) 

_ 1 f°° ^(^ + - 

^ ■'0 {v + Ai)2(t; + M2 + uyJ{v + M2 + uy - Auv 



(70) 
(71) 



and Af is given in Eq. (58). This integral can be evaluated and the result is given in Eq. (|A3|) 
of the appendix. 

The result in Eq. (^4]) and the more detailed study of Ref. [0, which obtained 
= 0.091 GeV, suggest that meson-loop effects should not make a large contribution to 

the normalisation. This is indeed what we find, as illustrated in Fig. |^ which shows that for 

reasonable values of Mc and A2: 



0.2 GeV < Mr < 0.6 GeV and 0.3 GeV < Ao < 0.7 GeV 



(72) 



the correction is small and negative, leading to a less than 2% correction at the physical 
value of the pion mass. 

The relative importance of the quark core and pion-loop contributions to the electro- 
magnetic pion charge radius, Tj^, can also be obtained from Eq. 



Substituting we obtain 
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f2 

J TT 



du 



[4m + 



L^214 



9W 



(73) 



(74) 



where (r^)*^^"^ is the generalised impulse approximation contribution. 
The result in Eq. (|3^ 



and that of = 0.59 fm obtained in Ref. suggest that meson- 
loop effects can additively contribute 10-15% to the charge radius at the physical value of 
the pion mass. The same conclusion can be drawn from the phenomenological Coulomb 
gauge DSE studies of Ref. This is just what we find, as we illustrate in Fig. |^ which 
shows that with f.,^ = 0.093 GeV, reasonable values of and A2, Eq. (0), and at the 
physical value of the pion mass, the correction to the square of the charge radius is always 
positive and less than 15%. 

It is of interest to study the chiral limit. It will be recalled that we have set the mass of 
all the external pions to zero so that the contribution of the loop-pions in the chiral limit 
is easily identified: — >■ 0. As we have remarked, the generalised impulse approximation 
contribution, (r^^^), is regular in the chiral limit and only weakly dependent on the current 
quark mass. In Fig. |^ we plot the ratio f'^/{r^^^)'^ as a function of for = 0.22 GeV and 
A2 = 0.5 GeV. This figure emphasises that at ^0.14 GeV the dominant contribution 
to the pion charge radius is provided by the quark core. It is not until becomes very 
small, on the order of 0.01 GeV, that the pion cloud contribution becomes as important as 
the quark core and this contribution is well described by the form 



[r. 



2\GIA 



0.73 - 0.082 In 



[m. 



(75) 
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for ml < 0.14 GeV. 

It is natural to compare our result with that of ChPT. In Ref. 
the result 



one finds in Eq. (5.6) 



V2LI 



1 



327r2 F2 



2 In 



+ In 



+ 3 



(76) 



where /i^ is the loop regularisation scale and Lg is one of the ten standard low energy 
constants in the effective action of ChPT into which an infinity from the divergence of the 
pseudoscalar loop has been absorbed [^ . 

The ten low energy coefficients, of which L'g is one, are fixed by requiring that the effective 
action of ChPT provide a good description of low energy strong interaction phenomena 
while the loop regularisation scale //^ is usually taken to be something of the order of 
or because this is a measure of the mass scale at which additional excitations should be 
accounted for. Of course, the actual values of the coefficients L[ depend on /x^ because of 
the ambiguities associated with regularising the divergent loops but, once a scale is chosen, 
values can be quoted. 

Fitting K% decay using /z^ = [Ref. 0, below Eq. (14.1)] and Fq = 93.3 MeV a value 

of 



TV 

^9 



(7.1 ± 0.4) X 10" 



(77) 



is quoted in Ref. in which case 

12 a 



(0.38 ± 0.020) fm^ 



(78) 



which is 90% of the experimental value: 0.439 ± 0.008 fm [Q. A similar calculation in 



Ref. PSi using /i = ml with 



LI = (6.9 ± 0.2) X 10"^ 
yields [upon correction of a minor typographical error therein] 



12L^ 
Fi 



(0.37 ± 0.011) fm^ 



(79) 



(80) 



which is 84% of the experimental value. One sees that at the physical value of the pion mass, 
and with an accepted value of the renormalisation scale, the "chiral logarithm" contributes 
little to the charge radius of the pion. 

The pion charge radius has also been studied in nonlinear chiral theories, Ref. |l|], where 
a comparison is made between the contribution of nucleon and pion loops. In this study 
the pion loop contribution at the physical value of the pion mass was directly identified and 
found to be 0.06 fm^ which is just 14% of the experimental value. 
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VI. SUMMARY AND CONCLUSIONS 



Herein we have calculated the charge radius of the pion in the coupled Dyson-Schwinger- 
Bethe-Salpeter equation approach as formalised in the Global Colour-symmetry model, 
which is a QCD based model field theory. In this approach, the "bare hadrons" are non- 
pointlike objects with a "quark core" obtained as the solution of Bethe-Salpeter (mesons) 
or relativistic Faddeev (baryons) equations whose kernel is chosen in such a way as to model 
the quark-quark interaction in QCD. The meson cloud associated with a given hadron is a 
small correction in this framework and is obtained through a systematic expansion based on 
the effective action of the GCM. 

In this approach the charge radius of the pion receives a contribution from its quark core 
and from pion loops. We showed that the quark core contribution is finite in the chiral limit 
and that, at the physical value of the pion mass, it is the dominant determining characteristic 
of the pion with the pion loop contribution being a small, finite, additive correction of less 
than 15%. The fact that the loop contribution is finite is a general property of this formalism; 
it is due to the internal quark core structure which provides a natural cutoff in all integrals 
that arise. Our result is consistent with a recent reanalysis of the pion charge radius obtained 



in lattice simulations 44 



We explicitly identified the origin of the chiral divergence in the charge radius of the pion 
as arising from the pion loop; this feature is not absent in the general DSE-BSE approach 
but is simply an higher order correction, as are all meson loop effects. Our calculation 
demonstrated that the pion's own pion cloud is not a very significant component unless the 
mass of the pion is less than ~ 50 MeV. 

A broad implication of our result is that it demonstrates that the many quark-based 
models which do not reproduce the logarithmic divergence of the charge radius in the chiral 
limit; for example, the relativistic, constituent quark model of Ref. |^5[, nevertheless capture 
the most important characteristic of the pion: its quark core, and should not be discarded 
on this basis alone. 

In closing we note that the contribution of the quark core of the proton/neutron to its 
charge radius has not been calculated in the approach we have used herein, although such a 
calculation is underway. In our opinion, and if the analysis of Ref. is a reasonable guide, 
the pion loop contribution will be a little more important in this case but, at the physical 
pion mass, will still generate less than a 30% correction. 
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APPENDIX: MISCELLANEOUS EXPRESSIONS 

In this appendix we present expressions for (rl)'^^^, Eq. and g{u), Eq. (|A3D, since 

they are rather lengthy. 

The generalised impulse approximation contribution to the pion charge radius is [O: 



where 



H{x = q^) = B [-A^B^B' K" K^x^ - A^B" K' K^x' - 36 A^B'K' 

-A^B' K" K^x^ + A^B A' A" K^x^ -30 A B^A' 

-42 A'^B'^B' - 24 A^B^K" K^x - 48 A'^B K' K^x 

-54 A B^A' K' K'^x - 42 A^B'^B' K' K^x - 72 A^B A' K^x 

-18 B^A''^K^x - 24 A^B' K^x - A2 A B'^A' B' K^x 

-24 A B'^A" K^x - 21 A^B^B" K^x - 4 ^^53^/3^2 

+4 A^B^K" K'Kx^-6A B^A' K'^K x^ - 2 A^B^B' K''^K x^ 

-24 A'^B K" K^x^ -3 A B^A' K" K'^x^ - 68 A^B A' K' K'^x^ 

-2 B^A'^K' K^x^ - 28 A'B' K' K^x^ - 18 A B^A' B' K' K^x^ 

-8 A B^A" K' K'^x^ + 3 A^B'^B" K' ifV - 20 A^B A'^K^x"^ 

-40 A^A' B' K^x^ -UABA' B'^K^x"^ + 2 A'^B'^K^x'^ 

-36 A^B A" K^x^ - 5 A" K^x"^ -AAB^B' A" K^x^ 

-8 A'^B" K^x^ + 2 A B'^A' B" K^x^ + 4 A^B B' B" K^x^ 

-4 A^B K'^x^ + 4 A^B K" K' K x^ ~ Q A^B A' K'^K x^ 

-2 A^B' K'^K - 3 A^B A' K" K^x^ - 10 A^B A''^K' K'^x^ 

-10 A^A' B' K' V - 4 A^B A" K' K^x^ -QAB A'^K^x^ 

-6 A'^A'^B' K^x^ - 2 A^B' A" K^x^ - 2 A^A' B" K^x^) . (A2) 



As we remarked above, the integral in Eq. ( [7lD can be evaluated. Writing x = M^/ K\ 
and M = u/Aj one obtains: 

[x — 1)^ 

^^^^ ~ ((1 + x) Ina; + 2(1 - - 1 + xf[{u - I + xf + Au] ^ ^^^^ 
[u + 1 — x){u — 1 + x)^ 
^Aux{u -1 + x) + (m+1- x){u + l + x)[{u-l + xf + Au] 



'u - 1 + a;)2 + 4u 
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X In 



[u - 1 + xY + - {u + 1 - x) 



(u + x)^J {u — 1 + xy + Au+ (u + x)^ + U — X ^ 

[(u + l)(u -l + x) + {u + l- x){u + l + x)][{u - 1 + xf + Au\ 



X In 



2\lu[u + x) 



^u{u + x) - 
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FIGURES 



FIG. 1. This figure is a pictorial representation of tlie amplitude identified with the tttt^ vertex 
in a generalised impulse approximation. The straight external lines represent the incoming and 
outgoing TT, the filled circles at the vr legs represent the {'n\qq) Bethe-Salpeter amplitudes, the 
wiggly line represents the photon, 7, the shaded circle at the 7 leg represents the regular part of 
the dressed quark-photon vertex [which satisfies the Ward-Takahashi Identity, Eq. (|^)] and the 
broken internal lines represent the dressed quark propagator. 

FIG. 2. This figure is a pictorial representation of our approximation to the pion loop contribu- 
tion to the 7r-7r7 vertex, Eq. (p^). The straight external lines represent the external pions, the large 
filled region at the tt legs represents T^^^\p + q,p — k,k — q), Eq. (^), the internal lines represent 
the virtual loop pions, the smaller filled region represents A^^^{k,p -\- q — k), Fig. and the wiggly 
line represents the photon. 

FIG. 3. In this figure we plot [/tt — fn], in MeV, in order to illustrate the dependence of the 
pion-loop induced renormalisation of /tt on Mc and A2, physically reasonable ranges for which are 
0.2 GeV < Mc < 0.6 GeV and 0.3 GeV < A2 < 0.7 GeV. 

FIG. 4. In this figure we plot [f'^/{r^)^^^ — 1] in order to illustrate the dependence of the pion 
loop contribution to (r^) on Mc and A2, physical ranges for which are 0.2 GeV < Mc < 0.6 GeV 
and 0.3 GeV < A2 < 0.7 GeV. 

FIG. 5. In this figure we plot [f'^ / {r"^)^^^ — 1] as a function on m^, in MeV, in order to illustrate 
the onset of the pion-loop induced In 111^^ divergence of (r^). It is clear from this figure that even 
for m-TT = 50 MeV the pion-loop contributes < 25% of the charge radius. The dashed line is the fit 
of Eq. (|7|). 
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